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Abstract 

Regularity and uniqueness of weak solution of the compressible isentropic Navier- 
Stokes equations is proven for small time in dimension N = 2,3 under periodic 
boundary conditions. In this paper, the initial density is not required to have a 
positive lower bound and the pressure law is assumed to satisfy a condition that 
reduces to 7 > 1 when N — 2,3 and P{p) — ap 1 . In a second part we prove a 
condition of blow-up in slightly subcritical initial data when p £ L°°. We finish by 
proving that weak solutions in T N turn out to be smooth as long as the density 
remains bounded in L 00 (L < - N+1+£ ^) with e > arbitrary small. 

1 Introduction 

The Navier-Stokes equations are the basic model describing the evolution of a viscous 
compressible gas, As emphasized in many papers related to compressible fluid dynamics 
|4 H H5 1 IBTfl 1621 163] . vacuum is a major difficulty when trying to prove global existence and 
strong regularity results. As a matter of fact, starting from initial densities that have 
positive lower bounds, local existence of smooth solutions can be proved by classical 
means, since lower bounds on the density persists for small enough time. This paper is 
devoted to the proof of well-posedness results. We want to prove next that the norm 
L v {L q ) on the pressure P(p) control the breakdown of strong solutions of the Navier- 
Stokes equations when N = 2,3. In other words, if a solution of of the Navier-Stokes 
equations is initially suitably smooth and loses it regularity at some later time, then the 
maximum norm of the density grows without bounds at the critical time approaches. 
Let us first recall the periodic compressible isentropic Navier-Stokes equations in 
(N > 2). 

(d t p + div(pu)=0 mP'((0,T)xT*), 

\ d t (pu) + div(/m ® u) - div(2p(p)Du) - V(A(p)divn) + VP(p) = pg, 

The unknowns p, u respectively correspond to the density of the gas p > and its 
velocity u £ M. N . The last equation of defines the pressure P which is assumed to 
be increasing in p for physical reasons. Usually P(p) = Pj >a = «P 7 for some positive 
constant a and some 7 > 1. The viscosity coefficients ares assumed to satisfy p > 0, 
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iVA + 2/z > 0, and the external forces g to belong to L 2 ((0,T) x T^)^ for all T > 0. 
Finally we complement the above system with initial conditions 

fp*=o = Po>0, 
[ pu t =o = m . 

Before the remarkable work of P-L Lions, very little was known about solutions of the 
compressible isentropic Navier-Stokes equation at least when N > 2. In [50] he proved 
a global existence theorem and weak stability results for Py a pressure laws under the 
following assumptions on the initial data 



( p eL\T N )nW(T N ),p >o, 



1 G L\Y N ), 
Po 



(1.3) 



where we agree that — ^ = on {x G T such that po(x) = 0}. More precisely he proved 



Theorem 1.1 We assume JO)) and 7 > 1 if N = 2, 7 > § if N = 3. Then there 
exists a solution (p,u) G L°°(0, 00; L 7 (T Ar )) x L 2 (0, 00; -ff 1 (T Ar )) Ar satisfying in addition 
p G C([0,c»),ZP(T")) i/1 <p < 7, p|n| 2 G L°°(0, 00; L 1 (T iV )), p G L<^[0, oc); L^)) 
for 1<(7<7 — 1 + ^7. Moreover, when f = 0, for almost all t > 0, we /iai>e 



-p\u\ 2 -\ —p" / )(t,x)dx + / (is / OulVuP + (A + /x)(divu) 2 )dxds 

2 7-1 Jo 



< / {^Po\u \ 2 + —^—Po)(x)dx. 
Jjn 2 7-I 



(1.4) 



In addition, he proved similar results for more general pressure laws P(p) such that 

- 1 P(s) 



, ds < +00, 

S (1.5) 
P(s) 

liminf -> s ^ +Q o — — > 0, 

for some 7 satisfying the above condition of theorem 11.11 Notice that the main difficulty 
for proving Lions' theorem consists in exhibiting strong compactness properties of the 
density p in L^ oc (R + x W N ) spaces required to pass to the limit in the pressure term 
P(p) = ap" 1 . Let us mention that Feireisl and his collaborators in [28, 29, 30J generalized 
the result to any 7 > ^ f° r N > 2 in establishing that we can obtain renormalized 
solution without imposing that p G Lf oc (R + x R N ) (a property that was needed in Lions' 
approach in dimension N = 2,3 giving the further condition 7 — 1 + ^ > 2), for this he 
introduces the concept of oscillation defect measure evaluating the loss of compactness. 
In [7J Bresch and Desjardins show a result of global existence of weak solution for the 
non isothermal Navier-Stokes system assuming density dependence of p(p) and X(p), 
considering perfect gas law with some cold pressure close to the vacuum, and the following 
relation: 

\(p) = 2(pp'(p)-p(p)). (1.6) 
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The key point in this paper is to show that the structure of the diffusion term provides 
some regularity for the density thanks to a new mathematical entropy inequality. This 
one has been discovered in [8], we call it the BD entropy. Mellet and Vasseur by using 
the BD entropy, get in [52] a very interesting stability result. The interest of this result 
is to consider conditions where the viscosity coefficients vanish on the vacuum set. It 
includes the case p(p) = p, X(p) = (when N = 2 and 7 = 2, where we recover the Saint- 
Venant model for Shallow water). The key to the proof is a new energy inequality on 
the velocity and a gain of integrability, which allows to pass to the limit. Unfortunately, 
the construction of approximate solutions satisfying: energy estimates, BD mathemati- 
cal entropy and Mellet- Vasseur estimates is far from being proven except in dimension 
one or with symmetry assumptions, see |54j . |49| . [32J. Note that approximate solutions 
construction process has been proposed in [5J satisfying energy estimates and BD math- 
ematical entropy. This means that only global existence of weak solutions with some 
extra terms or cold pressure exists in dimension greater than 2. 

The existence and uniqueness of local classical solutions for (jl.ip with smooth initial data 
such that the density po is bounded and bounded away from zero (i.e., < p < po < M) 
has been stated by Nash in [55]. Let us emphasize that no stability condition was re- 
quired there. 

On the other hand, for small smooth perturbations of a stable equilibrium with constant 
positive density, global well-posedness has been proved in [51] . Many works in the case 
of the one dimension have been devoted to the qualitative behavior of solutions for large 
time (see for example |4(J|, 148] ) . Refined functional analysis has been used for the last 
decades, ranging from Sobolev, Besov, Lorentz and Triebel spaces to describe the regu- 
larity and long time behavior of solutions to the compressible model [UT], [63J, [42], [45] , 

S3- 

The use of critical functional frameworks led to several new well-posedness results for 
compressible fluids. In addition to have a norm invariant by (jl.ip . appropriate func- 
tional space for solving (jl.ip must provide a control on the L°° norm of the density (in 
order to avoid vacuum and loss of ellipticity). For that reason, we restricted our study 
to the case where the initial data (po,«o) and external force / are such that, for some 
positive constant p: 

JV — — 1 — — 1 

(po -p)€ n G B;i, and / G L} oc (R + , G ) 

with (p,pi) G [1, +oo[ good choose. The most important result come from R. Danchin in 
[22] which show the existence of global solution and uniqueness with initial data close from 
a equilibrium, and he obtains a similar result in finite time. The interest is that he works 
in critical Besov space ( critical in the sense of the scaling of the equation) . More precisely 

JV — — 1 — — 1 

to speak roughly, he get strong solution with initial data in B 2 2 1 D B 2 2 1 x (B 2 2 1 ) N ■ 
Here compared with the result on Navier-Stokes incompressible, he needs to control the 
vacuum and the norm L°° of the density in the goal to use the parabolicity of the mo- 
mentum equation and to have some properties of multiplier spaces. That's why Danchin 
works in Besov spaces with a third index r = 1 for the density, and it's the same for the 
velocity as the equations are linked. In [25], R. Danchin generalize the previous result 
with large initial data on the density. 

In [25], however, we hand to have p = Pi, indeed in this article there exists a very strong 
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coupling between the pressure and the velocity. To be more precise, the pressure term is 
considered as a term of rest for the elliptic operator in the momentum equation of 
This paper improve the results of R. Danchin in [22^ I25j . in the sense that the initial 

JV 

density belongs to larger spaces i3 p p 1 with p G [1, +00 [. The main idea of this paper is 
to introduce a new variable than the velocity in the goal to kill the relation of coupling 
between the velocity and the density. In [TT], F. Charve and R. Danchin and in [13] 
Q. Chen et al generalize the results from [23J by choosing more general initial data. In 
particular they works with general Besov space constructed on L p , however they added 
some conditions on p {p < 2N) to get global solutions. For results of strong solutions 
with general viscosity coefficients we refer to [151 E31 G2] ■ 

In |36| . we address the question of local well-posedness in the critical functional frame- 
work under the assumption that the initial density belongs to critical Besov space with 
a index of integrability different of this of the velocity. We adapt the spirit of the re- 
sults of [1] and [33] which treat the case of Navier-Stokes incompressible with dependent 
density (at the difference than in these works the velocity and the density are naturally 
decoupled) . The main idea of this paper is to introduce a new variable than the velocity 
in the goal to "kill" the coupling between the velocity and the density. We introduce a 
new variable v\ to control the velocity where to avoid the coupling between the density 
and the velocity, we analyze by a new way the pressure term (in particular we will use 
this variable v\ in this article). This idea is inspired from the works of D. Hoff, P-L 
Lions and D. Serre about the the famous effective pressure. More precisely we write the 
gradient of the pressure as a Laplacian of the variable v\, and we introduce this term in 
the linear part of the momentum equation. We have then a control on v\ which can write 
roughly as u — QP(p) where Q is a pseudodifferential operator of order —1. By this way, 
we have canceled the coupling between v\ and the density, we next verify easily that we 
have a control Lipschitz of the gradient of u (it is crucial to estimate the density by the 
transport equation). This result allows us to reach some critical initial data in the sense 
that we are not very far to choose (po — Pi ^0) m 1 x v 

On the other hand there have been few existence results on the strong solutions for the 
general case of nonnegative initial densities. The first result was proved by R. Salvi and 
I. Straskraba. They showed in [59] that if f£ is a bounded domain, P = P(-) G C 2 [0,oo), 
po G H 2 , uq G Hq n H 2 and the compatibility condition: 

1 

Lu + VP(po) = pig, for some g G L 2 , (1.7) 

is satisfied, then there exists a unique local strong solution (p, u) to the initial boundary 
value problem (jl.ip . H. J. Choe and H. Kim proved in [17J a similar existence result 
when Q is either a bounded domain or the whole space, P(p) = ap 1 (a > 0, 7 > 1), 
po G L 1 n H 1 n W 1 ' 6 , uq G -Dq n D 2 and the condition (fOD is satisfied. 
B. Desjardins in [26j proved the local existence of a weak solution solution (p, u) with 
a bounded nonnegative density to the periodic boundary value problem (|1.1|) as long as 
SUP <KT*(I|P(*)IU°°(T3) + l|Vu(i)|| L 2 (X 3)) < +00. 

This paper is devoted to improve the works in |26] and [IT] by choosing very low regularity 
on the velocity. In the sequel we will note 4 = dt + u • V and f = 4+f ■ 
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The viscosity coefficients will suppose constant in the sequel and are assumed to satisfy: 

5 

p>0, 0<A-/i. (1.8) 
It follow that there is a I > 4, which will be fixed throughout, such that: 

^>i^ (19) 
X > 4(1-1)' (y) 

In the sequel we will assume that g G E\ with: 

\\g\U = \\g\\iw) + WqWlw + f U{s) 7 \\v f\\ 2 L , + f f f( s ) 5 \d s f\ 2 dsdx, 

T T T Jo Jo Jr n 

where f(s) = min(l,s) . We obtain now the following existence of weak solutions in 
finite: 

Theorem 1.2 Let N = 2,3. Assume that p, and A verify il.8\) and \1.9\) . We assume 
i i 

that po G L°°(R N ), PqUq G L 4 and PqUq G L 2 . Moreover g is in E\, and g G L 4 (L 4 ) n 
L°°(L 4 ). 

• There exists Tq G (0, +oo] and a weak solution (p, u) to the system Al.l}) in [0, To] 
such that for all T < Tq and with f(t) = min(i, 1): 

sup / ^p{t,x)\u{t,x)\ 2 + \P(p(t,x))\ + f(t)\Vu(t,x)\ 2 dx 



0<t<T 



+ sup / [lp(t,x)f(t) N (p\u(t,x)\ 2 + \Vu(t,x)\ 2 )dx (l.io) 
o<t<T Jrn I 

r+oo n 

+ / [\Vu\ 2 + f{s)p\u\ 2 + \uj\ 2 ) + a N \\7u\ 2 }dxdt<CC . 

Jo Jr n 

sup / p(t,x)\u{t,x)\ 4 dx < C . (l.ll) 
o<t<T Jm N 

where Co depends of the initial data po and uq. We obtain moreover: 

pe L°°((0,T) xR N )nLf(Ll(R N )), (1.12) 



(1.13) 



{ ^d t uEL 2 (L 2 (R N )), 
VtVu G L 2 T (H 2 (R N )), 

VtG = Vi[(X + 2/x)divu - P(p)} G Lt(H 1 (R n )), 
^VuG Lf{L 2 (R N )), 
where V denotes the projection on the space of divergence- free vector fields, and 
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In addition if we assume that uq G H 2 1+e with e > 0, we obtain the following 
estimates: 

sup / [~p(t, x)\u(t, x)\ 2 + \P(p(t, x))\ + t 2 "^ - e \Vu(t, x)\ 2 dx 
*>o Jm^ ^ 

+ sup / [~p(t,x)t a (p\u(t,x)\ 2 + \Vu(t,x)\ 2 )dx (1.14) 

0<t<T Jr^ ^ 

+ / / [|Vu| 2 + i 2 --->| 2 + ^Vu| 2 ]dxd£<C(Co + C/) e , 
Jo Jr n 



where: 



N 

a = 3---e, ifN = 2, 

N 3 
cr = max(3-y-e,6--iV-3e), if N = 3, 



(1.15) 



sup HV) || ^_ 1+e <CC« and sup ||^«(t, < C(r)C 9 , (1.16) 

Kt<T 

u>ii/i r G (1,4] and: 



0<t<T f?7r " 0<t<T 



Vti G L^(BMO). (1.17) 

• TTie regularity properties 111.10)) , kl.lS\) . {1.1$ , 111.16}) and \l.lty hold as long as: 

SU P \\p\\l?°(L°°(RX)) < +°°- ( L18 ) 

*6[0,T] 

Remark 1 • In this theorem we do not assume hypothesis on the vacuum as in 
\3J$ where — G L°° . However we get a control ^(BMO) on the gradient of the 
velocity. We know that for incompressible Navier-Stokes this hypothesis is enough 
to get uniqueness. In this sense we can consider our result as a theorem of strong 
solutions in a weak sense. 

• In this result we improve the works of B. Desjardins in [26] by the fact that we 
can choose more general initial data. The second important point is that estimates 
111. lb]) and Ill.lH) can hold as long as p G L°° in dimension 2 and 3. 

Let p > 0. In the sequel we will note q = p — p. In the following corollary we obtain 
strong solution by adding some regularity on the initial density and we assume that po 
is away from the vacuum. The goal is to get by this extra regularity on the density a 
control Lipschitz of Vu. 

Corollary 1 Under the hypothesis of theorem \1.2L we assume in addition that there 
exists c > such that po > c. Moreover we assume that qo G B^ ^ if P{p) = Kp and 
qo G #2Voc if P is a general pressure. The solutions of theorem \1.2\ are then unique and 
verify locally in time \1.10\) . i!.13\) . fi.i^p , i!.16\) and: 

VuG L^ (L°°). 

Moreover if p in L°°(M. x T N ) then the solution is global. 
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In the following theorem we want improve the criterion of blow-up of corollary [TJ More 
precisely we prove that it is just necessary to control the norm L°°(L( N+1+e >i) with e > 
of the density when P(p) = ap' 1 with 7 > 1 to obtain global strong solution. This result 
is to connect with the works of Serrin for incompressible Navier-Stokes system where in 
the compressible Navier-Stokes, the pressure plays the role of the velocity. 

Theorem 1.3 Let A = and g as in theorem HOI and g G L°°{L°°). Let P(p) = ap 1 
with a > and 7 > 1. Assume that (p ,u ) G (L^ n L°° n x (I? n L°° n H^~ 1+t ) 

with e > 0. Moreover po check po > c > 0. 

Let (p, u) a weak solution of system il.l]) on [0, T) with the previous initial data which 
satisfies the following condition: 

p€ L oo (0,+oo,L( iV+1+ ^), (1.19) 

with e > 0. Then (p, u) is unique and verify locally in time tl.lO\) . M.13]) . {1.1$ , M-lOjl 
and: 

Vu g L} oc (L°°). 

Remark 2 • Ln this theorem we can see that by compare with the incompressible 
Navier-Stokes equation, the good variable to control is not the velocity but the pres- 
sure. Indeed if we control enough the pressure, we get integrability on the velocity. 

• In this theorem we need to assume hypothesis on the vacuum as in \3J$ where 

G L°° in the goal to get a control Lipschitz on the velocity. 

• This initial data are considered slightly surcritical in dimension 3 except that uq G 
L°° (it is probably possible to improve this fact). 

• Here we need to assume that A = to get a control L°° on u as in the article of A. 
Mellet and A. Vasseur in [53]. We recall that in this paper they need of a control on 
P{p) G L°°(L 3+£ ) with e > for N = 3 by using some De Giorgi technics used by 
A. Vasseur in t° reprove the famous result of Caffarelli-Kohn and Nirenberg in 
F77J/. As in 153]/ , the pressure plays a important role, and in some sense the pressure 
is the good variable to control to get global strong solution. The pressure plays the 
role of the velocity for Navier-Stokes incompressible when we have compressible 
Navier-Stokes system where a structure of type effective pressure exists. 

• We can justify the previous statement by seeing the results of D. Bresch and B. Des- 
jardins in fdjj. They have a new mathematical entropy for a class of such viscosity 
coefficients which gives some norms on the gradient of p (see |6j/) ; in particular 
we can obtain the relation U.19\) . However this type of viscosity coefficient kill 
the structure of effective pressure and we can apply our proof. It's again a proof 
that the structure of the viscosity coefficients plays a crucial role for compressible 
Navier-Stokes system. 

• As in the case of incompressible Navier-Stokes system (see ]10$). there is a big gap 
between obtaining A1.19\) and by this way have a control on p in L°° . 
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• V. A. Waigant has builded in [66] explicit solutions for which the maximal integra- 
bility of the density correspond to L g (0,l, L q ) with q = 7 ^ 3JV 2 ^^ — , It means that 
h!.19\) fails in this case except that the force term is less regular than in our case. 
It means that the regularity of g is crucial to get strong solution. 

Remark 3 We believe that our method can be adapted to the euclidian space R . This 
is the object of our future work. 

Our paper is structured as follows. In section [21 we give a few notation and briefly 
introduce the basic Fourier analysis techniques needed to prove our result. In section [3] 
and HI we prove a priori estimate on the density and the velocity. In section [5] and section 
[6l we prove the theorem 11.21 and corollary [TJ We finish in the section [7] by the proof of 
theorem 11.31 



2 Littlewood-Paley theory and Besov spaces 

Throughout the paper, C stands for a constant whose exact meaning depends on the 
context. The notation A < B means that A < CB. For all Banach space X, we 
denote by C([0, T],X) the set of continuous functions on [0,T] with values in X. For 
p G [1, +oo], the notation L p (0, T, X) or Lj,(X) stands for the set of measurable functions 
on (0,T) with values in X such that t — > \\f(t)\\x belongs to L p (0,T). Littlewood-Paley 
decomposition corresponds to a dyadic decomposition of the space in Fourier variables. 
We can use for instance any ip G C°°(T N ), supported in C = G T N < |£| < |} such 
that: 

5>(2- J o = i if e/o. 

Denoting h = J r ~ 1 ip, we then define the dyadic blocks by: 

Aiu = (p(2- l D)u = 2 lN / h(2 l y)u{x - y)dy and Siu = V A k u . 

Formally, one can write that: 

u = ^ A k u . 

This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in 5 (T ) for two reasons: 

1. The right hand-side does not necessarily converge in S' (T N ). 

2. Even if it does, the equality is not always true in S' (T N ) (consider the case of the 
polynomials). 

2.1 Homogeneous Besov spaces and first properties 
Definition 2.1 For s € K, p G [1, +oo], q G [1, +oo], and u G ^'(T^) we set: 

Hk g = (£( 2, l A Hk*) 9 )'- 

The Besov space Bp q is the set of temperate distribution u such that \\u\\b» < +oo. 
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Remark 4 The above definition is a natural generalization of the nonhomogeneous Sobolev 
and Holder spaces: one can show that B^ ^ is the nonhomogeneous Holder space C s and 
that .Bf 2 is the nonhomogeneous space H s . 



Proposition 2.1 The following properties holds: 

1. there exists a constant universal C such that: 
C~ 1 \\u\\b s < II Villi < C||ii||_B3 . 

ii ll-"p,r n "-Dp,r 11 " p,r 

2- If Pi < P2 and n < r 2 then B s pi Ti J Bp 2 "^ {1/pi_lM) . 



3. Bp n <—> Bp r if s > s or if s = s and r\ < 



r. 



Before going further into the paraproduct for Besov spaces, let us state an important 
proposition. 

Proposition 2.2 Let s£R and 1 < p, r < +oo. Let (u q ) q >-i be a sequence of functions 
such that 

(^2<HKHl P )"<+oc. 

<7>-l 

If suppui C C(0, 2 q Ri, 2 q R 2 ) for some < R\ < R2 then u = Ylq>-i u g belongs to Bp r 
and there exists a universal constant C such that: 

IHI^ r <c 1 +H(^(2^|K|| LP r)i 

<?>-! 



Let now recall a few product laws in Besov spaces coming directly from the paradifferen- 
tial calculus of J-M. Bony (see [1]) and rewrite on a generalized form in [1] by H. Abidi 
and M. Paicu (in this article the results are written in the case of homogeneous spaces 
but it can easily generalize for the nonhomogeneous Besov spaces). 

Proposition 2.3 We have the following laws of product: 

• For all s£l, (p, r) £ [l,+oo] 2 we have: 

\\uv\\ B b t < C(\\u\\ L °°\\v\\Bi r + \\v\\L°°\\u\\B* ir ) ■ (2.20) 

• Let (p,pi,p2,r,\i,\ 2 ) e [l,+oo] 2 such that:± < ^ + ±, p 1 < \ 2 , p 2 < M, ~ < 
— + and - < i + j^. We have then the following inequalities: 

if 8! + 82 + NM(0, l- 2 ±-±)>0, s 1 +» < K anda3+ » < M. then: 

\\uv\\ _ n( m + j s .x ) < h\\ B ;i Jv\\ B s2 i0o , (2.21) 
when Sl + g = ^ (resp s 2 + f = we replace WuWb^JHb^ (resp IH^J 

h W u \\B;iJ v \\B a P i r ( res P IHbz^l^)' *M + £ = |r and s ^ + f = » we take 
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r = 1. 

// Sl + s 2 = 0, si G (f - — , - - y- ] and i + i < 1 t/ien: 

' 1 V Ai P2 ' Pi A 2 J Pl P2 — 



|ra|| .^(i + i.i, < |M|| B n ||w|| B »2 . (2.22) 

o n Pl + P2 P J " B P" " 



'PI P2 P' PI i 1 P2' 

|s| < ~ for p > 2 and — 4- < s < 4r e/se, we /iaue: 
p p p ' 



|uu||rs < CIIuIIrs IMI a . (2.23) 



Remark 5 In the sequel p will be either p\ or p2 and in this case -3- = — — if p\ < P2, 



Corollary 2 Let r £ [1, +oo], 1 < p < p\ < +oo and s such that: 
• -e(-^+JV(l + i-l) > g)^I + i>l J 

JV 

i/ien we /iaue if u £ Bp r and v G -Bp^oo n -L°° : 

||w||b s < C|l u ll-B s IMI 



A pi P2 



p,r — ' " p,r ■■ ■■ pi 



RP1 nroo 

-'-'p^jOO 1 ^ 



The study of non stationary PDE's requires space of type L p (0,T, X) for appropriate 
Banach spaces X. In our case, we expect X to be a Besov space, so that it is natural to 
localize the equation through Littlewood-Payley decomposition. But, in doing so, we ob- 
tain bounds in spaces which are not type L p (0, T, X) (except if r = p). We are now going 
to define the spaces of Chemin-Lerner in which we will work, which are a refinement of the 
spaces L^Bp^). 

Definition 2.2 Let p G [l,+oo], T G [l,+oo] and si G E. We set: 

ii u ii^(B-) = (E 2Zrsi ii A ^)ir^(L P) )"- 

zez 

W^e then define the space L P n (B^ r ) as the set of temperate distribution u over (0,T) x 
such that \\u\\r P m s 1 •> < +oo. 

We set CT(Bp^ r ) = L^(Bp l r ) t~) C([0,T], B^ r ). Let us emphasize that, according to 
Minkowski inequality, we have: 



Remark 6 It is easy to generalize proposition \2.3l to Lj,{Bp\) spaces. The indices s\, 
p, r behave just as in the stationary case whereas the time exponent p behaves according 
to Holder inequality. 
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In the sequel we will need of composition lemma in Lj,{B pr ) spaces. 
Lemma 1 Let s > 0, (p,r) G [l,+oo] and u G L^(B^ r ) n L^(L°°). 

1. Let F G W i 2 +2 '°°(T JV ) suc/i too* F(0) = 0. Then F{u) G L p T (B s pr ). More precisely 
there exists a function C depending only on s, p, r, N and F such that: 

ll F ( U )IIZ£(B^.) ^ C ^\ U \\L^(L^))\W\\~ L P T{B s r y 

2. Let F G W^ ] C +3 '°°(T N ) such that F(0) = 0. Then F{u) - F'(Q)u G L p T {B s pr ). More 
precisely there exists a function C depending only on s, p, r, N and F such that: 

\\F(u)-F'(0)u\\ lPT{B s r) < C(IMlL~(L~)||u||fp (B| , r) . 

Here we recall a result of interpolation which explains the link of the space B p 1 with the 
space -Bp j00 , see [2T] . 

Proposition 2.4 There exists a constant C such that for all s G R, e > and 1 < p < 

+oo, 

« \\ZP( B b )<C \\U \\lP (Bs Jl + log-jpj: . 



Now we give some result on the behavior of the Besov spaces via some pseudodifferential 
operator (see [2T]). 

Definition 2.3 £e£ m G R. A smooth function function f : T N — > R is said to be a S m 
multiplier if for all muti-index a, there exists a constant C a such that: 

V£eT N , 1^/(01 <c a (i + ieir-w. 



Proposition 2.5 Lei m G R and f be a S m multiplier. Then for all s G R and 1 < 

p,r < +oo i/te operator f(D) is continuous from B pr to B p ~ m . 

We now focus on the mass equation associated to (jl.ip 

I d t q + u- Vq + qh(q) = -qdiw 1} 
\ Q/t=o = Qo- 



N 

where h G C°°, h(0) = and // G W S '°°(R,R). Here «i belongs in L 1 {B p p ^) with e > 
and p G [1, +oo]. 
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Proposition 2.6 Let 1 < p < p\ < +00, P2 G [1, +00] and 1 < r < +00, 1 < r\ < +00. 
Lei assume that: 

11 JV 
-JVminf— , — ) < <r < — , (2.25) 
Pi P P2 

with strict inequality if r < +00. Assume that q G 5£ r , Vui G L\Q,T; L°°), divt>i G 

JV 



Lj,(Bpi i0 a) H L°° and i/ia£ g G L 1 ^ {B° r ) satisfies (2.24)- There exists a constant C de 



pending only on N such that for all t G [0, T] and m G Z, u>e /zave; 

ll9llzr(^)^ eW( °ll*lkr + ^ w -i. (^6) 

with: 

ft AT 

V(t)= {\\Vu(t)\\ m +[|divt;a(r)|| jv + h^Wlt 1 ) dr if a <l + -, 

Jo Bp^coni™ B/.f i00 nz°° Pi 

= / (IIVu(t)H jv + ||div«i(r)|| jv + \\q(r)\\l+ l )dT if a < 1 + - 

Jo Bg >00 nL<x> Bp P 2 2 oonL°° Pi 

and r = 1. 

wnt/i a i/ie smallest integer such that a > s. 
Proof: Applying A; to ([272^]) yields: 

d t Aiq + u • VAiq + Ajg = Ri - Aj(gdiv«i) - Ai{qh(q)) with flj = [it • V, Aj]g. 
Multiplyng by Aiq\Aiq\ p ~ 2 and performing a time integration, we easily get: 

\\Aiq(t)\\ LP dT < \\Aiq Q \\ LP + / (H-R/JIlp + ||divn||Loo||A^|| LP 

Jo 

+ ||Ai(gdivi;i)||iP + ||A z (g%))|| LP )dr. 
By paraproduct, there exists a constant C and a positive sequence (q) G T such that: 
||Aj(gdivui)[|Lf. < Cci2~ la \\q\\ Bz ||divui| 



R p 2 nroo 



Similarly by lemma [T] we have: 

||a,(«%))||£p ^Cc^ikll^JI^II^ 1 - 

Next the term H-RjHlp may be bounded according to the inequality (2.53) p 110 of [3]: 

||(2 ZCT ||^|| iP ),|| r <C||Vn|| jv \\q\\B° r . 

l(—+e) 

We end up with multiplying the previous inequality by 2 p ; and summing up on Z: 

lk(*)b pV < lko||B pV + J*cv'\\ q \\ B « r dT + CV'dr. 

Gronwall lemma yields inequality (|2.26p . □ 
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3 A priori bounds on the density 

In this section we make a formal analysis on the partial differential equations of 
and begin by classical energy estimates. Multiplying the equation of conservation of 
momentum by u, we obtain 



1 



(-p\u\ z (t,x)+U(p)(t,x)dx + / (pD(u) : D(u)(s,x) 

JN 2 ./|| JjN 

+ (A + p)\divu\ 2 (s,x))dsdx < [ fi^-(x) + Il(p )(x)})dx, 



(3.27) 



jn 2p 



where II is defined by 



U(s) = s[ I ^-dz), (3.28) 



o 



It follows classicaly that we have the following bounds 

'2\ 



p £ L°°(0,oo;L^ 

L°°(0,oo;L 2 ), (3.29) 
VneL 2 (0,oo;Ly 2 . 



3.1 Bound on logp 



Let us emphasize that one of the main ingredient of the proof of the first part of theorem 
11.21 is a partial differential equation derived from (jl.ip involving log p. It was introduced 
by P-L Lions in [50] to prove global existence of weak solutions of (jl.ip in a particular 
case and it is one of the key of the proof in the paper of B. Desjardins in [26]. Letting 
formally (A) _1 div operate on the equation of conservation of momentum, we obtain 

{2p + X)divu - P(p) + P(p)dx = dtA" 1 div{pu) + R i R j (pu i u j ) (3.30) 

where A -1 denotes the inverse Laplacian with zero mean value on T N and Ri the usual 
Riesz transform. Let us observe that the equation of mass yields 

dt log p + n.V log p + divn = (3.31) 

Let us define F and G by the following expression: 

F = (2p J + X){logp + A- 1 div(pu)), 
G= (2p + X)divu- P(p). 

Moreover we shall denote respectively by P and Q the projection on the space of 
divergence- free and curl- free vector fields. Combining (|3.30p and (|2.24p . we obtain 

d t F + u ■ VF + P(p) - [ P(p)dx = [uj,RiRj](pUi). (3.32) 

Next we define the Lagrangian flow X of u by 

(d t X(t,s,x) = u(t,X(t,s,x)), 
\ X /t=s = x, 
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and derive the following identity 

F(t,X(t,0,x)) = F (x)- [ P(p(s,-))dxds+ [ {[uj, RiRj}(pui)(s, X{s,0,x))ds. (3.34) 

Jo Jo 

Using the fact that p(-) > 0, we obtain 

F(t,x) < F (X(0,t,x))+ [ P(p(s,-))dxds + I {[uj,RiRj\{pUi)(s,X(s,t,x))ds. (3.35) 



It follows that 

\og{p{t,x)) <log(||po||L-) + C||(A)- 1 divm || L ^ +C||(A)- 1 div(H(t,-)||L- 

ft (3.36) 

+ C t + C / \\[uj,R i Rj](pu i )(s,-)\\L<x>ds, 
Jo 

where Co depends of the initial data. In view of the usual Sobolev embedding inequalities, 
we obtain 



log(p(t,a;)) <log(||po||L-) + C||(A)- 1 divmo||L-+C||(A)- 1 div(p«)|| Lo 



+ C t + C I {{[u^RiRj^fmi^s,-)]^! ds, 



(3.37) 



with e > 0. Let us now remark that we have 

||Vn|| L iv +e < C(||curl«|| £ w +e + \\G\\ l n+, + \\P(p)\\ L »+<)- (3-38) 
In view of R. Coifman, P-.L. Lions and S. Semmes [19], the following map 

W l ' ri (T N ) N x L T2 (Y N ) -> W 1,rz (T N ) N 



(a,b) -4 [aj,RiRj]bi 



(3.39) 



is continuous for any N > 2 as soon as ^ = ^j- + — . Hence we have the following 
estimates for N = 3: 

\\[uj,RiRj](flUi)(8, ^Wli^^ < C\\Vu\\ L 2 iL6) \\pu\\ L 2 iL 6+, ) , (3.40) 
with e > 0. We obtain finally: 

Iog(p(t,s)) <Iog(||po|U«) + C||(A)- 1 divm ||L- +C||(A)- 1 div(pu)|| L - 

+ C Q t + C\\Vu\\ L 2 (Lli) \\pu\\ L 2 {L&+t) , 

4 A priori estimates for the velocity 
4.1 Gain of integrability of the velocity u 

We want here derive estimate of integrability on the velocity u. This idea has been 
successively used in different papers, we refer in particularly to [33] and [52j[53]. To do 
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it, we multiply the momentum equation by u|u| Pl 2 and we get after integration by part: 
— f p\u\ Pl (t,x)dx + p f I (\ur- 2 \Vu\ 2 {t,x) + : ^^\u\ Pl - 4 \V\u\ 2 \ 2 (t,x))dxdt 

Pi JjN J JjN 4 

+ A / / ((divu) 2 \u\ p ^ 2 (t,x) + ^^divuy" Ui di\u\ 2 \u\ Pl ~ 4 (t,x))dtdx 
Jo Jjn 2 / — J 

~ [ I P{p){divu\u\ P1 - 2 + (pi - 2)V UiU k diU k \u\ Pl - A ){t,x)dtdx 
Jo JtN i, k 

< [ p \u \ p 'dx. 

We have then by Young's inequality: 

A(Pl - 2) f I divuV Uidi\u\ 2 \u\ pi - 4 (t,x))dtdx 
2 Jo Jt n i 

= \P±^[ [ divuu- V{\u\ 2 )\u\ pi - 4 {t,x)dtdx < 
2 Jo Jt n 

A P 1 _-2 ( r? /•* /■ | divu |2| u | P1 -2 dtd2 , + 2 /"* /" | V | n |2|2| u | P1 -4 ( ^ x) ^ dx) 
2 2 Jo Jt^ V Jo Jt n 

If we choose: 

for some s £ (0, 77), by the fact that (divii) 2 < iV|Vu| 2 we therefore obtain: 

— / p\u\ Pl {t,x)dx + A s I I \u\ Pl - 2 \Vu\ 2 {t,x)dtdx 
Pi Jt n Jo Jt n 

+ B S [ [ \u\ Pl ~ 4 \V\u\ 2 \ 2 (t,x)dxdt < [ [ P(p)(divu\u\ P1 - 2 
Jo Jt n Jo Jt n 

'"' ~ 2 -u ■V{\u\ 2 )\u\ p ^ 4 ){t,x)dtdx + f p \u \ Pl dx, 

JjN 



+ 



2 

Pl -2„ (pi-2) 2 A 2 



with A s = p(l — sN) and B(s) = ^^—p — i6(s^+A) • By Young inequality, we get again: 

- / p\u\ Pl (t,x)dx + A s f f \u\ Pl - 2 \Vu\ 2 (t,x)dtdx 
n Jt n Jo Jt n 

[ [ \u\ pl ~ 4 \V\u\ 2 \ 2 (t,x)dxdt <C e [ [ P(p) 2 \u\ pi - 2 dtdx+ I p \u \ pi dx, 
Jo Jt n Jo Jt n Jt n 



Pi 

+ B S 



with C e enough big. Now we want use the fact that Vjul^ 1 £ L 2 (L 2 ), which implies that 
when N = 3 u G L P1 (L 3pi ). More precisely we have: 

IMI^i (L 3 P1) < C(l|V(|n|^)|| L 2 (L2) +upx), 
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where upi) is the average of \u\ 2 . We have then by Holder's inequalities with p ^ p 2 + 
H%±11 = 1 and 2i=2 + JL = l: 

3pi pi pi 

P{pf\u\ pi - 2 dtdx\ < \\P(p) 2 \\ gi 3P1 |||n| pi - 2 || pi 3pi , 

Jj N L t ^"(L 2 (Pi+ 1 )) l pi- 2 ( L pT^) 

< ll^(^)ll 2 3 P1 IMIgT^y 

<c||P(p)|| 2 3M (||v(H^) + H^(-)|| i?(i2) ) 2 ^. 

Remarking that J^ N podx = M / 0, we can write as 7 > |: 

Pi 1 Ex f Pi_ 

\u\-(s)<—(\\p{s,-)\\ L ->\\V\u\^{s,-)\\ L 2+ / p(s,x)\u(s,x)\^dx), 

M JfN 
* ^H~ 2 I 1 pj_ 1 J_ £l 

u| 2 (s)ds) 2 < j^[\\p\\L^(Ly)\Mu\ 2 || L 2 (L 2) + {Mt)2\\pnu\\l r{LVl) ), 
We obtain then: 

II^(p)II 2 3P1 (l|v(H l2L ) + |^(-)|| i?(L2 )) 2 ^< 







3 P1 (I"'" 



Ci||P(p)|| 2 a Pl ([|V(|ti|*)[|.,? f « + (Mt)— ||^n|' P1 - 2 



By a standard application of Young inequality ( 2 ^~ + = 1), we obtain that: 

— / p|u| pi (t,x)dx + A, / |u| pl - 2 |Vn| 2 (t,x)dtdx 
Pi Jt n Jo 

+ B S f\ur- i \V\u\ 2 \\t,x)dxdt<C 2 et \\P{p)\\ 2pi 3P1 +- / 
where C £j t is big enough and depend of the time t 



(4.42) 

po\u \ Pl dx, 



4.2 Gain of derivatives on the velocity u 

In this section we deal with the case N = 3. The case N = 2 follows the same lines. In 
the sequel we will follow the procedure developped in and [35] to get some energy 
inequalities. The main idea compared with the results in [161 El HB] is to obtain energy 
inequalities which depends only of the control on p £ L°°. It implies that we have to 
be carreful to not introduce some derivatives on the density which means to "kill" the 
coupling between velocity and pressure. Multiplying first the equation of conservation of 
momentum by f(t)dtu with f(t) = min(l, t) and integrating over (0, T) x T N , we deduce 
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that: 

ft 



f f f(t)p\d t u\ 2 dxds + l [ f{t)(p\Vu(t,x)\ 2 + {X + p)\dwu\ 2 (t,x))dx 

JO JT N - JT n 

[ [ VP(p)- f(t)d t udxds< [ [ ^-(p\Vu(t,x)\ 2 + (X + p)\divu\ 2 (t,x))dxds 
Jo Jt n Jo Jt n 2 



rt 

+ 



+ / IIv7(*)p^IIl2(t^)(II\/p(^- V)u\\l»(tn) + \\Vpf\\v>crf))ds. 
Jo 

(4.43) 

Next we use the equation of mass conservation to write: 



t 

o Jt n 



VP(p) ■ f(t)d t udxds = -dt [ f(t)P(p)divudx + f d t (f(t)P(p))divudx 

JjN JjN 

= -d t [ f(t)P(p)divudx - [ f(t) [div(P(p)u)divu + (pP' (p) - P(p))divu] dx 

JjN JjN 

- / / (t)P(p)divudx, 

JjN 

= -d t [ f(t)P(p)dWudx + —^— ! f(t)P(p)u-V(G + P(p))dx 

JjN -/' • A JjN 

~ (2/ , + A) 2 J JN f^PP'iP) ~ P (P)){G 2 - P{ P f + 2(A + 2p)P(p)divu)dx 

- / f\t)P(p)divudx, 
Jtn 

If we define Hf(s) = s(j s ^-dz), we have then by mass equation: 

dtU f (p) + dib(n / (p)n) + /(p)divw = 0. 
By this fact we obtain that: 

P{p){u-VP{p)-2{pP\p)-P{p))divu) = P{p){-d t {P{p)) - 3P' (p)pdwu + 2P(p)dwu). 
and: 

f VP(p) • f{t)d t udx = -d t [ f(t)P(p)dwudx + -r^-dt [ f(t)k(p)dx 
+ I f(t)P( P )u • VGdx + - 1 — 9 I f(t)P 2 (p)( P P'(p) - P{p))dx 



2p 



T\ /„ '«> p M" • vafc + (vW /r» /W^W^W - p 

(2?w jC / ( ( ) g2 (" p '(") - p ^ - xh L '' mp)dx 

f'(t)P(p)divudx. 
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where k(s) = P(s) 2 — P(p) 2 + s(J^ —g^-ds—^^-). Inserting the above inequality in (|4.43p 
and by Young's inequality, we obtain: 

f f(t)p\d t u\ 2 dxds + \ [ f(t)(p\Vu(t,x)\ 2 + (A + p)(dwu(t,x)) 2 )dx 

JjN 2 JjN 



+-^— / f(t)k(p(t,x))dx<C+ [ f(t)P(p(t,x))dwu(t,x)dx 
A + 2p Jjn Jjn (4.44) 

+-!—[*[ f(t)k(p)dx + I I f'(t)P(p)dwudx 
A + Z V Jo Jt n Jo Jt n 

+C f f f(t)(\ P P'(p) - P(p)\G 2 + \P{p)u\\VG\ + \^pu- Vu\ 2 
Jo Jt n 

+\Vpf\ 2 )d*ds. 



In the sequel we set: 



A(t)= [ I f{t)p\d t u\ 2 dxds + \ [ f(t)(p\Vu(t,x)\ 2 + (X + p)(divu(t,x)) 2 )dx 

Jo JjN 2 JjN 

+ (0l 1 X s 2 f I f(s)P(p) 2 (pP'(p)-P(p))dxds + —l—[ f(t)k(p(t,x))dx 
(2p + X) z J J T N X + 2pJ T N 

We obtain finally: 

A(t)<C + C t \\p\\ L «, + C I [ f(s){\ P P' (p) - P{ P )\G 2 + \P(p)u\\VG\ 

Jo Jt n 

+ \Jpu ■ Vu\ 2 + \^pf\ 2 )dxds. 

<c+c Aii/ p'(p)-p( /0 )iiL-iiy7MGiii 2 +/( S )ii 5 (p)( S ,-)iiL-iiv / p^iiL 2 , A AK , 

Jo (4.45) 
x ||VG|| i2 + \\^pu\\l4y/mVuf L4 + \\p\\ L ™\\vW)f\\h)dxds, 

<C + C [\\\h(p(s, \\y/f(sjVu\\ 2 L2 + f(s)\\i(p(s, + f(s)\\VG\\ L i 



o 



x \\g(p(s,-))\\L°o + IIVP«lli*llv7WV«|li4 + HpllL-llvTM/lliO^, 

where g(p) = /i(s) = \sP'(s) - P(s)\ and i(s) = h(s)P(s) 2 . 
Estimates on Pu and G 

We want now to obtains bounds on Vu and g, assuming that p is a priori bounded in 
L°°(T ). Indeed we wnat show that the control of -A(t) in (|4.45p depend only of a control 
on HpH^oo. 

Next we use once more the equation of conservation of momentum to write: 

pAu = V(pd t u) + V(pu ■ Vu) - V(pf), (4.46) 
VG = Q(pd t u) + Q(pu ■ Vu) - Q(pf). (4.47) 
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Therefore we have: 

||VG|| L2 + \\AVu\\ L 2 < C\\p(s, Olllco (|| sfpd s u{s-)\\v + \\Jpu- Vu(s, -)\\ L 2 

The case N = 3 

For simplicity, we will treat only the case of the dimension 3. We recall that for all 
1 < p < +oo: 

||Vu||£* < ||Wu||iP + \\RG\\ LP + \\R(P(p))\\ LP . 
We want now to recall the Gagliardo-Nirenberg's theorem: 

V/a 1 ^) such that J N fdx = 0, ll/II^T^^CII/lll.^llV/llJ,^. 

We deduce that from Gagliardo-Nirenberg's inequality and Young's inequalities: 



[IVHll4[|V7WV«[|!« < Cf(s)\\^u\\U\\R(P(p))\\% + \\Wu\\ L * + \\RG\\ L ,) 
< C||^||i 4 (/( S )||P(p)||| 4 + (x/7(i)(||Vn|| L4 + \\P(p)\\ifi))t 



x(^m(\\AVu\\ L2 + \\VG\\ L i))>) 



< a(/( s )||^«||| 4 ||p( P )||| 4 + \\ P (s r )\\l^u\\l4^m(\\vu\\ L , + \\p(p)\\ l i))i 



< C(/( S )||Vpn||| 4 ||P(p)||| 4 + -f(s)(\\Vu\\l 2 + ||P(p)||! 2 )l|p(^-)llioo||^lli 4 



x (Vm(\\Vpdsn(s-)\\ L2 + \\^pu-Vu(s,-)\\ L2 + \\p(s,-)\\l ac \\f(s,-)\\ L2 )*) 

1 

+ ef(s)(\\^pdMs-)\\h + \\Vpu- v«(-, -)\\h + OIIl- OlliO) • 

(4.49) 

Hence we obtain by Young inequality: 

f(s)\\g( P (s, oili- I|vg( s , < f -)||i- ||p(*. Oil WOO (4 50) 

+ e(\\ ^fWpd s u(s-)\\h + ||>/tfto« • v«0», + / WIIp(-, Oil*- ll/0», Oil!*)- 

By adding (|4.49p and (|4.5Q[) . we obtain: 



IIVH&II y/m^\\h + /(«)ll5(p(s, 0IU- l[VGf( s , -)\\^ < 
c(f(s)\\^-pu\\l4P(p)\\ 2 Li + i/00(||v«|& + ||P(p)||i»)||p(*,-)lli-ll^lli* 
+ e /( s )(|| Vp^( s -)II| 2 + Olli— ll/(«, OHIO) + -/WIIM*. Olli- 

x l|p(a,-)l|i»- 

Therefore we have by using inequality (|4.42|) . (|4.45p and (|4,51|) : 

A(t) < c + c /%(IIp(OIIW/(*)((i + ||V< 2 + \\P(pM)\\h) 

+ WfM\\h)dB, 
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(4.51) 



o 



(4.52) 



where C depends of the time t here. Gronwall's lemma provides the following bound: 

A(t) < e X p(Cexp( /V(IIpOv)||l~ )<&)), (4.53) 
J o 

where <f> G C°(R+,1R+) n C 1 (0, oo) such that </>(s) > eo-s for some positive s. 
Control of sup 0<t<T f 2 (t) f \ii\ 2 (t, x)dx + / / f 2 (s)\Vu\ 2 dxds 

In the sequel, we want obtain estimate on Vu in L\{BMO), that's why we need of addi- 
tional regularity estimates. We derive then estimates for the terms f(t) 2 J JN \u\ 2 (t, x)dx 
and Jq J T N fN(s)\Vu\ 2 dxds. First we rewrite the momentum equation on the following 



r* 

form: 



pii - pAu - (A + p)Vdivu + VP(p) = pf. 

We apply to the momentum equation the operator ^ = dt + « • V, we recall here the 
following equalities: 

^pu j = p^h? + d t pvP + u J ^2 d k p u k , 

k 

= p—vP — pdivuvP , 
dt 



We have next: 



u—Au j = ud t Au j + ^ d k Au j u k , 

at -J 



dt 

k 



= e(dtAv? + div(Au J 'u) — AnMivu), 

and (where Z) = divu): 

(A + p)—djdivu = (A + p)(dtdjD + div(5jDu) — djDdivu), 
We obtain finally: 

P ju j + a,a t P(p) + diva,P(/>) = p(d t Au j + div(A^u) 

+ (A + p){d t djD + div(a jJ Du)). 
We shall make use of the following transport theorem if pw = f\ and if h = g(t), then: 

J I -d s (hpw 2 )dxds = I I (—h pw 2 + hwf)dxds. 
Jo Jt n 2 Jo Jt n 2 

We apply the previous result with: 

/i = -djdtP(j>) - dwdjP(p) + p(d t Au j + dw(Au j u) + (A + p){d t djD + div(0j£>u)), 
and with /i(s) = f(s) 2 , we obtain then: 

' • " 2 [ n (+ <r\U,(+ rA\ 2 ^ — [ [ — -p ( c\ f'( a\ nXiifi rlnrrl a -L f [ f(e\ 2 ;j 

(4.55) 



-f(t) 2 p(t,x)\u(t,x)\ 2 dx= / -/(«)/ ( S ) / o|«| 2 dxd S + / / /(s)^ 
2 J Jo Jtn 2 J J t n 



x [-(djd t P + dw(djPu) + piAdtu 1 + div(A^u)] + (A + p)[djd t D 

+ div((9j-Z)-u)](ix<i,s. 
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Since f(s)f'(s) < f(s) we can apply (|4.53p to bound the first term on the right. Next 
by integrations by part we get: 

-// f{s) 2 u j (d j d t P + div{d j Pu))dsdx= f f f(s) 2 (d j u j d t P + d k u j d j Pu k dsdx, 

Jo J7 N " ' •/[) JjN 

= J J f(s) N P' \dju j d t p + d k u j dj(m k )dsdx 



o Jt n 



f(s) 2 P'[-d j u j (pd k u k + d kP u k ) + d k u j d jP u k }dxds, 



s,x)\ dxds. 



= - I [ f{s) 2 [P' pDdjU j + d k Pu k djU j - djPu k d k u j }dxds, 
Jo Jt n 

= ~ f f f{s) 2 [P pDdjU j - P(DdjU j - dju k d k u j )]dxds. 
Jo Jt n 

This term is therefore bounded in absolute value by: 

c( f I p{ P )f{s) 2 \vu\ 2 dxd S y{f [ f{ s ) 2 \vu\ 2 dxd S y 

Jo Jt n Jo Jt n 

<C e [ [ P(p)f(s) 2 \Vu\ 2 (s,x)dxds + e [ [ f(s) N \Vii( 
Jo Jt n Jo Jt n 

The third term on the right side of (|4.55p may be written: 

-M / [ f 2 [Vit j -Vul + iVrf -^Au^dxds 
Jo Jt n 

= -li\ [ f 2 [Vu j ■ (Vi4 + V(W ■ u)) + u\(u k vj u - (uju l ) k )]dxds, 
Jo Jt n 

= -p / / 2 [Vii j • (Vu{ + V(Vu j ■ it)) + u j k (u k u 3 u - uj k u l - uju l k )]dxds, 
Jo Jt n 

2|V7 ,V 1 2 



<-/W / f\Vu 3 r + M f 2 \Vu\ 2 (\Vii\ + \D\)dxds. 

JO JjN J J 

The last term on the right side of (I4.55P may be bound as follows: 

~{X + p) f I f 2 {djd t D + dw{djDu)dxds <-(A + p) I f f 2 D 2 dxds 
Jo Jt n Jo Jt n 

+ M [ [ f 2 \Vu\ 2 (\Vii\ + \D\)dxds. 



JJ 



A" 



It then follows by Young's inequalities that: 



f(t) 2 [ p\u\ 2 (t,x)dx+ [ [ f 2 (s){p\X7u\ 2 + (X + p)\D\ 2 )dxds 

JjN ./|| JjN 

ft 



<M[C + C £ Co\\P(p)\\l°°+ [ [ lf(s)f'(s)p\i 

Jo Jjn 2 



u\ 2 dxds 



+ 11 f(s)\Vu\ A dxds], 
Jo Jt n 



(4.56) 



<M[C + C e Co\\P(p)\\ L °°+A(t)+ f [ f 

Jo Jt n 



2 (s)\Vu\Uxds]. 
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Next from the momentum equation, we obtain as in the works of D. Hoff in []: 
p\Vw\ 2 = p(div(wVw) + dj(pwit k ) — dk(pwii^) + p{v? dkW — u k djiu). 
Integrating we thus obtain: 

/ I f(s)\Vw\ 2 dxds < M( I I f(s)p\u\ 2 (s,x)dsdx+ [ [ f(s)p\w\ 2 {s,x)dsdx), 
Jo Jjn Jo JjN Jq JjN 



< M{ f [ 
Jo Jr N 



f(s)p\ii\ (s,x)dsdx + 



L° 



<M(A(t) + \\p\\ L ~). 

(4.57) 

Similarly we obtain: 



sup f(t) 2 \Vuj\ 2 dx < M( / f 2 (t)p\u\ 2 (t,x)dsdx + f(t) 2 p\w\ 2 (t,x)dx) 

0<t<T J JjN JjN 

<M([ f 2 (t)p\u\ 2 (t,x)dsdx + C\\p\\ L °o), (4.58) 

JjN 

<M(C + C e C \\P(p)\\ L o B +A(t)+ f [ f 2 (s)\Vu\ 4 dxds + C\\p\\ L oo). 

Jo Jt n 

To complete the estimates (|4,56p and ()4.58p . we will need of estimation on JjJ J TN f 2 (s)\Vu\ (s, x) A dxdt 
and to conclude we will use (|4,53p . 

Control of f Q f TN f 2 (s)\Vu\(s, x) 4 dxdt 
We have then: 

[ [ f 2 Wu\ 4 dsdx< I f f 2 ({RG) A +u A ){s,x) + f 2 {s)RP{pf{s,x))dxds. (4.59) 
Jo Jt n Jo Jt n 

Let focus us on the case N = 3. When iV = 3we can apply Gagliardo-Nirenberg, let: 

pG|| L 4<C||G||| 2 ||VG||| 6 . 
We have then by Young's inequality and the fact that (A + 2p)AG = div(p(u + g)): 

f f f 2 ( s )(RG)' l ( S ,x)dsdx< f f 2 ( S )\\G\\l 2 \\VG\\l G ds, 
Jo Jt n Jo 



o 



* 5 

L 2 1 



< / 2 (^)l|G||^(||Vu|| i2 + ||/|| L6 )^ S , 



t r ft 

2/' t ,MY7„\|2 i n / &( a \\\m\\W . 



<e/ / f 2 (s)\Vu\ 2 dxds + C e / f 2 (s)\\G\\%ds, (4.60) 
Jo Jt n Jo 



<ef f f 2 \Vu\ 2 dxds + C e [ (f( s )\\G\\v) 10 -±^d 8 , 

Jo Jt n Jo f{ s ) 

<e/ / f 2 \Vu\ 2 dxds + C e C( sup /(s)||G|| L2 ) 10 . 

70 Jt n 0<s<t 
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From (|4.60p and (|4.56|> . (|4.58|> . we can conclude that: 

B(t)f(t) 2 [ p\ii\ 2 (t,x)dx + f I f(s)(p:\Vu\ 2 + (X + f,)\D\ 2 )dxds 

JjN J JjN (4.61 J 

< M[C + C e C ||P(p)|U» + CC e A(t)(l + A(t) 9 )}. 

sup f(t) 2 [ \Vuj\ 2 dx<M(C Q + C t C \\P(p)\\ L oo + A(t)C\\p\\ L oo 
o<t<T J (4.62) 

+ eB{t) + CC e A{t) 10 ). 

We can remark that all the inequalities (|4.53|) . f)4.57j) . (|4.6ip and (|4.62p depend on the 
control of ||p||i<x>. In the following subsection, we want explain that we can control p in 
L°° in finite time. 

Conclusion 

We will treat by simplicity only the case N = 3. We want here to explain how to control 
the norm L°° of the density p in finite time. From ()3.3T[) . we have: 

\og(p(t,x)) <log(|| / oo||L-) + C||(A)- 1 divm ||L- +C||(A)~ 1 div( / on)llL- 

+ C \\[uj,RiRj](pUi)(s,-)\\ L oods, 
Jo 

From the previous section and the regularizing effects on u we obtain: 

log(p(t,x)) <C t + Cexp(Cexp(/ <j>(\\p(s, )<**))■ 

J o 

Consequently by Gronwall lemma, there exists Tq > such that for all T < Tq: 

\\p\\l°° < C. 

Proof of (|Q3|) and ([TTB|) 

We want get now solutions such that we control Vu in L {BMO). To do it, we need of 
additional regularity on the velocity. We will use again the technics introduced by D. Hoff 
in [39]. The idea is to obtain some estimates by interpolation in "killing" the coupling 
between pressure and velocity. First , we mollify initial data satisfying the conditions 
of theorem 11.21 and then appeal to the result of [36] to otain a solution (p, u) defined at 
least for small time. We want now derive estimates on the solution independent of the 
mollifier and depending only on the initial data. 

In the sequel we will treat only by simplicity the case N = 3. Fixing the local in time 
solution (p,u) described above on the interval [0,To] with To > 0, we therefore assume 
throughout this section that C^ 1 < p < C with C > 0. We define a differential operator 
£ acting on functions w : [0, T] x —> T N by 

Cw = dt(pw) + div (pu ® w) — pAw — AVdivu>, 
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and we define w\ and w 2 by: 

Cw 1 = 0, (w 1 )/ t=0 = u , Cw 2 = -VP(p), (io 2 )/t=o = 0. (4.63) 

We observe here that by uniqueness w\ + w 2 = u. Straightforward energy estimates then 
show that: 

sup / \w 1 (t,x)\ 2 dx+ [ [ \Vw!\ 2 dxdt <C [ \u \ 2 dx (4.64) 

and: 

f-To 

\w2(t,x)\ 2 dx + / 

0<t<T Q Jl N 



sup / \w 2 (t,x)\ 2 dx + [ ° f \Vw 2 \ 2 dxdt < CT sup \P(p(t,-)\ 2 . (4.65) 
ct<T Jt" Jo Jt n o<t 



We shall derive (|1.14p and (|1.16p simultaneously as consequences of estimates for the 
following quantities: 



0<t<T 

for k = 0, 1 and, 



sup t^ k [ \Vwi{t,x)\ 2 dx + / / t 

<t<T JjN J J T JV 



1 k *wi\ 2 dxdt, 



sup / \V w 2 {t , x)\ 2 dx + i i 
o<t<i 7 Jo J 



w 2 \ 2 dxdt. 



To derive these bounds, we multiply equation (|4.63p for w\ and w 2 by w\ and 102, re- 
spectively and integrate. The details which are nearly identical to those in the previous 
section are quite straightforward, the essential point being that the spatial gradient of p 
must be avoided. Indeed the procedure of D. Hoff in [39J allows to "kill" the coupling 
between the velocity and the pressure. The results are that with C > 0: 

\p{r k f \Vw 1 (T,x)\ 2 dx\y T = t + f f T k p\wi\ 2 dxdr< 

J Jo Jr» (4 g6) 



1 

2' 



I T k - 1 \Vw 1 \ 2 dxdr+ [ 
o Jt n Jo Jt n 



^-pk I / T k ~ l \S7wi\ 2 dxdT + / / T k \Vwi\ 2 \Vu\dxdr, 



and 

ft 



-p f | Vw 2 (t, x)\ 2 dx + I j p\w 2 \ 2 dxdr < 

2 Jt n Jo Jt n ^ qj\ 

P(p(t,-))divw 2 (-,x)dx\l+ / (\X7w 2 \ 2 \Vu\ + \X7w 2 \\X7u\)dxdT, 

Jo Jt n 



By proceding exactly as in the previous section, we obtain the following results: 
/ \Vwi{t,x)\ 2 dx + / 

0<t<T , 



sup / \Vw 1 (t,x)\ 2 dx+ [ [ \wi\ 2 dxdt < C\\u \\ 2 H i, (4.68) 
\<t<r Jt n Jo Jt n 



sup t[ \Vwi(t,x)\ 2 dx+ [ [ t\wi\ 2 dxdt < C\\u \\ 2 L 2, (4.69) 
o<*<t Jt n Jo J 
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sup / \Vw 2 (t,x)\ 2 dx+ [ [ \w 2 \ 2 dxdt <CC , (4.70) 
o<t<T Jt n Jo Jt n 

where Cq depends only of the initial data. Now since the solution operator uq — > w\(t, ■) 
is linear, we can apply a standard Riesz-Thorin interpolation argument to deduce from 
(USED and IjOgjt that: 

sup / \Vw 1 {t,x)\ 2 dx + / / t x -P\wi\ 2 dxdt < C||u ||^. (4.71) 

0<t<T J Jo Jj N 

As u = w\ + w 2 , we then conclude from (|4.7U|) and (|4.7ip that: 

sup t 1 -^ I \Vu(t,x)\ 2 dx+ [ ° [ t l -P\u\ 2 dxdt<CC \\uo\\ 2 H p. (4.72) 
o<t<T Jt n Jo Jt n 

The next step is to obtain bounds for the terms 

sup i 2_/3 f \Vu(t,x)\ 2 dx+ ( I t 2 ~ p \Vii\ 2 dxdt 
o<t<T Jt n Jo Jt n 

appearing in (j!.14j) . To do this, we multiply the momentum equation of (jl.ip by t 2 ^^ii 
and integrate. The details are exactly as in the previous section, except now we apply 
the ft dependent smoothing rates established in (|4.7ip . Combining these bounds with 
(gjgj) , (1^65]) and K7T\i . we then obtain (fLT4"D for times t < T . To prove (fTTTB]) . we 
observe that for k = 0, 1, 

sup \\wi(t, -)\\ Hk < C\\u \\ H k, 

0<t<l 

by (ITOjl and (|Q5|I . Thus: 

sup \\wi(t, -)\\ H p < C||« || H fc, 

0<t<l 

for ft G [0,1]. As u = tfi + ^2, and applying ()4.70p we obtain that: 

sup \\wi(t, -)\\ H p < CC , 

0<t<l 

and then for r £ (2, g^j) in the case that ft > 0, that: 

sup ||ti(r, •) — < CCo- 

0<t<l 

This proves (|1.16|) . 

Regularity on the gradient of the velocity u in (W l,a + BMO) with a > iV 

Here we want examine the regularity of the gradient of the velocity and prove that Vu 
is in L\, q (BMO) to prove (|1.17p . More precisely we will verify that the new variable 
v\ introduced in [36] called "effective velocity' belongs in L)p g (W 2 ' a ) with a > N, let 
V^i G L 1 Tq(L°°). We recall here the definition of v\ introduced in [36]. The idea of |36| 
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was to introduce a variable v\ which allows to kill the coupling between the velocity and 
the pressure in the momentum equation of (jl.ip .In this goal, we need to integrate the 
pressure term in the study of the linearized equation of the momentum equation. To do 
this, we will try to express the gradient of the pressure as a Laplacian term, so we set: 

Av = VP(p). 

We have then v = (A) _1 VP(p) with (A) -1 the inverse Laplacian with zero value on T^. 
In the sequel we will set: 

X + 2u 
v i = u . 

V 

We have then: 

A«=VF + divw + (2^ + A)- 1 V(P( / o)), , x 

V ' (4.73) 

= Av 1 + (2 f i + X)- 1 V(P(p)). 

We can easily show that f^ \\Vv\\i°°dt < +00 if ()1.14[) holds. To see this, we apply 
standard elliptic theory on v\ and the fact that AG = div(pn — pg). We will use in 
particular the fact that divui = (X + 2p)G and curlt>i = w.We consider here only the case 
N = 3. The case N = 2 follows the same lines. For some a > 3 and e > determined 
by a we have then: 

IIVwxHl- < C(||VF|| L a + \\Voj\\ L c), 

< C(\\pu(t,-)\\ La + \\pg\\L" + ||Vw||i«), 

< C(\\p\\ L co\\u(t, •) - ii\\ L l \\Vu(t, + ||/3|U°°||fl'||i« + ||Vw||l<*), 

l-e e l+e 

<c(||p|U«»(||«(t,-)IIJ +u 2 )||v«(t, -)llxJ + IIpIMMU" + l|Vw|| L <»), 

We recall here that pii = 'pg, we have then as - p > C on [0, To]: 



u — Tr~n — P9- 

\\P\\ 1/2? 



so that by using (|1.14|) . we obtain: 

f ° \\Vvi\\ L ~dt < C [ V(i 1_s [ \u\ 2 dx)^r{t a f \Vu\ 2 dx)^dt + C , 
Jo Jo Jt n Jt n 

<C([ °&dt)*+C , 
Jo 

with s = y + e — 1 (e > 0) and where 4/3 = (s — 1)(1 — e) — (a + e). The above integral 
is therefore finite as 2/3 > —1. A similar result result holds for N = 2 with s > 0. Thus 
for the solution constructed in the previous section, ||Vi?i(i, OIU 00 *^ is finite if (|1.14p 

holds, inf p > c > and uq G P^ +€_1 ,with e > 0. 
More precisely we have proved in fact that: 

Vv x G L^ 1 '") Z^(i^). (4.74) 

with a = N + 2e where e > 0. 

As P{p) G L°° we deduce from (|4.73p and the results of C alder on- Zygmund, that: 

Vu G L\ (BMO). 



26 



5 Proof of theorem 11.2 



The above arguments are not rigorous, since we have to assume that p(t, x) is positive 
for all (t,x). In order to deal with possibly vanishing densities, we remark that if we 
assume as in [36J that po is also bounded from below, we can get L°° bounds for log p. 
In that case, when N = 2, 3, there exists Tq > such that for all t < Tq: 

Hp|ll,° o ((0,t)xT Jyr ) + II — H J L°°((0,t)xT iv ) - Ct- (5.75) 

Thus p is also bounded from below for small enough times, so that vacuum does not 
form on [0, To]. It follows that starting from a general bounded initial density p$, we can 
apply the above arguments to a weak solution (p n ,u n ) corresponding to initial values 
Po = Po + h an d Uq = uq which converge strongly in L°°(T N ) to po and uq. In view 
of the weak stability results of system (jl.ip given by E. Feireisl et al in [30], p n and 
y/fFu n respectively converge to p and u in C([0, T ],L q (T N )) and L 2 ((0,T ) x T^) for 
all q < 7 — 1 + Hence the uniform L°° bounds on p n yield L°° bounds on p. 
In the above formal derivations, we assumed that there exists global weak solutions of 
(jl.ip . This problem does not occur when we take 7 larger than y via the works of Feireisl 
et al in [3D]. When 1 < 7 < y, we can approximate solutions of (|l.ip by a global weak 
solutions of (jl.ip by a global weak solution (p n , u n ) corresponding to a modified pressure 
law that satisfies: 

P n (p) = P(p) + -p 2 , (5.76) 
n 

and the same initial data (pq,uq). Applying the above arguments on (p n ,u n ), we obtain 
all the uniform bounds for p n and u n on [0, To], where To does not depend on n.As a result, 
we also have uniform T 2 ((0, T ) x T N ) bounds on Vu n and L°°(0, T ,L 2 (T N )) bounds on 
yff^u n . Hence the weak stability results hold since the initial data are not functions of n 
and (p n ,u n ) converge to (p, u) in L 2 ((0,Td) x T n ) n+1 , where (p, u) is a solution of (jl.ip 
on (0,To). We refer to [50j for complete details of the stability proof. Let us emphasize 
that one of the key arguments for proving weak stability of solutions of (jl.ip is to obtain 
uniform L 2 ((0,To) x T ) bounds on p n to renormalize the transport equation. This is 
the case here as we have uniforms bounds on the density in L°°((0,To) x T^). 

6 Proof of corollary [L] 

6.0.1 Control of p G L°°(^ oc ) 

In view of proposition 12.61 where in our case h{p) = P(p), we have for all t € [0, T] and 
< e < 1: 

\\Ph r( B Uoo )< eCm ^+\\Po\\B^J, (6.77) 

where V{t) = f (||Vn(r)||ioo + [|diwi (t) H^g^ ^ + ^(r)!!^)^, where s the smallest 

integer such that P' £ W s >°°. We have seen by P~7lD that V«i G 1/(0, T, £^ i00 ) with 
e > 0, the main difficulty is to control Vu G 1/(0, T, L°°), for this we recall that: 

l|V«||Li ( r«») < l|V^i||^ (BSooo) + llPC^Ilii^). 

< II V«i h^B^) + WpWl^d^Wph^B^)- 
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Unsurprisingly the result comes from the well known following estimate: Next we have: 

\\p(t)\\ B o <C pt ) B o x Jog(e+ *V — , 
and we recall the following inequality: 

(5 1 
Vx > 0, V5 > 0, log(e + -) < log(e + -)(l + log<5). 

x x 

We obtain then from the previous inequality 

\\p(t)\\B0^ < ||p(t)|| B o ooo (l + log(||p(i)||^ oc ))log( e + 1 ^ 

Let X(t) = J Q \\p(s)\\ B o ^ds, we have then: 

t 



V(t)<C(X(t)+ / {\\Vn{r)\\ L ^ + \\^yvi(r)\\ BSB Jdr). 
Jo 

Combining (|6.77|) and the previous inequality leads to: 



X(t) < I \\p(s)\\ B<L (1 + CV(t) + log(l + lipoid J) Iog(e + rr rJ )ds, 

< / 11^)115^^(1 + CX(t) + C / (||V«i(r)|Uoo + 11^1^)1^ .Jdr 

JO j0 

+ log(l + UpoIIb^ oc )) log(e + t—t ttt: )ds. 



v 711 J -'oo,oo 



Applying Gronwall inequality and inequality (|4.74p shows that: 

X(t) < C tj0 exp(C^ \\p(s)\\ BOoa oa log(e + )cfe), 

<C 4 , exp(C /"*(!+ ||p(«)||£=o)dfl), 

JO 

where Ctfl depnds only of the time t and the initial data. As p G Lf^L 00 ), we conclude 
that X(t) < Ct and by this way we have proved that: 

IMU-OB^) < Co,t, (6.78) 
where C^o depnds only of the time T and the initial data. 

6.0.2 Control of p G ^(B^) 

In this case, we need to show for the sequel that p G L°°(Bj^ 1 ), and for this we proceed 
exactly as previous. Indeed as po G -B^r 1; by proceding as in the previous section we can 
show that p G L ( ^{B\ } 1 ). Next by using proposition 12.61 we obtain the fact that: 

llpll^^^e^Cl+lkollB^), (6.79) 

with V(T) = J T (||V«i(r)|| B i + e + [|ff(r)[|~ )dr. 
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6.1 Uniqueness 



We now discuss the uniqueness of the solutions of theorem ll.2I For this we want use the 
result of P. Germain [31] which is a result of weak-strong uniqueness. In the sequel we 
will note (pi,u\) the solution of the theorem 1 1 . 2 1 which exits on the time interval [0, To]. 
We have shown that our solution check p G L 00 (L 00 ). By theorem ll.2[ we obtain that 
our solution verify the following inequalities: 

sup / [-p(t,x)\u(t,x)\ 2 + \P(p(t,x))\ + a(t)\Vu(t,x)\ 2 dx 

0<t<+oo Jr^ 2 

+ sup / ip(t,x)f(t) N (fi\u(t,x)\ 2 + \Vu(t,x)\ 2 )dx 

o<t<+oo.M 2 (6.80) 

r+oo r 

+ / / [|Vu| 2 + f{s)p\ii\ 2 + M 2 ) + a N \Vu\ 2 ]dxdt 
Jo Jt n 

< C(C + C f ) 9 , 

and we obtain moreover: 



{ s/pdfU G L 2 (L 2 (T N )), 
VtVu G L 2 T (H 2 (T N )), 

s/tG = Vt[(X + 2p)divu - P(p)} G L 2 T (H l (Y N )), 
k ViVu£ Lf{L 2 {T N )), 



(6.81) 



Now by the result of P. Germain in [31] . we are able to prove that (p, u) = (px,u%) 
on [0,To]. To see this we have just to verify that verify the conditions of the 

theorem 2.2 of [31]. For simplicity we prove only the result for N = 3. We know that 
Vpi G L°°(B e Nl ) ^ Lf Q (L N ) and Vui G L 1 Tq {L°°). The main thing is to see that 
Vtui£L 2 To (L N ). 

We recall by Gagliardo-Nirenberg inequalities that: 

vt\\Ul\\ L 3 < 2 ||n 1 || i2 )2 (t2 a ][Viii||£a)2ta . 
From the inequalities (j!.14p . we deduce that \ftu\ G L 2 (L 3 ). 



7 Proof of theorem 11.3 

7.1 How to obtain a regularizing effect on v% when p E L°°(L g ) 

For the moment, we don't give conditions on q, we will precise his value in the sequel of 
the proof. 

We want now use our change of variable v i introduced in the previous sections. The 
interest of this new variable is to "kill" in a certain way the coupling velocity-pressure. 
In this goal, we can now rewrite the momentum equation of system (jl.ip . We obtain 
then the following equation where we have set v = 2p + A: 

pdtu + pu ■ Vu — /iA(-u v) — (A + /u)Vdiv(ti v) = pg, 
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where we recall that v = (A)" 1 (VP(p)) with (A) -1 the inverse Laplacian with zero mean 
value on T N . As v \ = u — we have: 

pdtV\ + pu ■ Vu — pAvi — (A + p)Vdivv\ = pg pdfV. 

v 

As d\vv = P(p) — Jj N P(p)dx, from the transport equation we obtain: 

divdtv = —P(p)pdivu — VP(p) • u + P' (p)pd\vu + VP{p) ■ u 

= -dxv{P(p)u) + (P(p) - pP'{p))divu - P(pjd\vu + pP'(p))divu. 

In the sequel we will need to use the Bogovskii operator that we note A -1 (see [SB] pl68 
for a definition), we obtain then: 

d t v = A- 1 (- div(P(»u) + (P(p) - pP'{p))dxvu - P(p)d\vu + pP'(p))divu). (7.82) 
We get finally 

pdtV\ — pAv\ — (A + ^)Vdivt>i = pg — pu ■ Vu pdtv. (7.83) 

We set f(t) = min(t, 1) and we remark that f(0) = 0. We multiply then (|7.83|) by 
f{t)dtV\ and integrating on (0,T) x we obtain where £ = p + A: 

f f{s)p\d s v 1 \ 2 dxds + \ I /{^(plVv^x)] 2 + ^diYv 1 ) 2 {t,x))dx < 

JjN I JjN 

f'(s)(p\V Vl (t, x)\ 2 + ^div^) 2 ^, x))dxds (7.84) 
+ / / pu ■ Vu f{s)dtv\dxds + / / {pg pd s v)f(s)d s vidxds. 

Jo JjN Jq JjN V 

We have next to control the terms on the right hand side of (|7.84p . In the sequel we 
will for simplicity treat only the case N = 3. We can recall now that from the works of 
Mellet and Vasseur in [53J , we control the velocity u in L°° (L°° ) as we have assume that 
p € L°°(L 3 ^ +e ) with e > 0. In fact from the inequality (|Q2|) , we could get a gain on 

1 3p 3jp 

ppu G L°°(L P ) with p arbitraly big if P(p) £ L p (Lp+ 1 ) and in this case p £ L 1p (Lp+ 1 ). 
We begin now with: 



t 

o Jt n 



t 

o Jj n 



u-Vupf(s)d s v 1 dxds= [ [ (u ■ Vv 1 + u ■ V(A) _1 V(P(p))) pf{s)d s v 1 dxds. 

Jo Jt n 



Estimates on the term f* J JN u ■ V(A) x V(F(p)) pf{s)d t Vidxds 
We start with treating the easiler term: 

\[ [ u-V(A)- 1 V(P{p))pf(s)d t v 1 dxds\< 
Jo Jf N 

Ce J J f(s)p\u ■ V(A)- 1 V(P( /9 ))| 2 dxd S + e\\^fWpdtVi\\ 2 L 2 iL2{TN)) 
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q_ 37 — q q 

As p E L°°(L q ) and so P(p) E L°°(i> ), we have p 1 u £ L°°(L 3 ^-^ ) (we recall that here 

q < 3 7 ). We have then to treat y/(i)p^^p~|n||V(A)- 1 V(P( / 9))| in (L 2 (T iV ). We 
must have in this case by Holder's inequalities: 

3g~67 7 37 - g 3g - 67 + 87g - 2g 2 1 2 , , 

2 +- + = 7T2 < o 6< 1 -2g 3 + 87) + 12 7 > 0. (7.85) 

2q z q q 2q z 2 

In this case, we get: 



t p pt 2q 39-67 

f(s)p\u-V(A)- 1 V(P(p))\ 2 dxds< / /( S )||p||2L7 Lg) || /0 || L ^)H' 9 III-(^) (is 



And we can conclude. 



Regularizing effect on Av\ 



We want here using the regularizing effect on v\ . To do it we use the momentum equation 
(|7,83p and we have: 



fiAvi + (A + //)Vdivui = pd t vi + pu ■ Vv 1 + pu ■ V(A) _1 V(P(p)) + pg - ^d t v. (7.86) 
We want use the ellipticity of (|7.86p to deduce regularizing effects on v 1 . 



Estimate on the term J Q || f(s)pu ■ Vv\\\ 2 L2 ds 



f(s)pu ■ Vvi\\ 2 L2 ds 







3q — 67 



< / Hp 9 -"II 2 9 \\Vpf( s )^ v i\\in \\p\\ Lq 9 ds, 







Here by Gagliardo-Nirenberg, we get: 

i^-Dw N#llT7 II 2 



1 e e_ 

gi 2 "T" 2q N 4q 2N 



with ^- = ^ + ^7 — 4f — x;— and q > By Holder's inequalities, we have to verify 



that: 

1 , J_ _ J_ _ i L + 3 ^~^ 1 3 g~ 6 ^ < j 

2 2? N Aq 2N q 2q 2 ~ ' 
It means that we must have: 

ll<7 2 -6<?(2 + 37) + 187 >0. (7.87) 
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We have next: 

\W f{s)pu ■ Vui \\ 2 L 2ds, 



< / ||p < u\\ 2 L ^\\Sns)Av x \\^ ~ 2) f{s)2\\Vv l \\ I l\\p\\L< l q ds, 







37-g 



< / ||p 9 «|| 2 9 (||p||2 9 ||v / p/('S)5twi|| L 2 + ^ff{s)\\pu-Vvi\\ L pds 

Jo 

+ ^m\\pu ■ v(A)- 1 v(P(p))|| LP + vW)\\pg\\L P + vW)\\pdtv\\ LP ) 2{1 -^ 

2- e 3q-6l 

x \Wvi\\jl {L2) f{s)2\\p\\ Lq q ds, 



< / 7>^P « /(s)l|V^|| i2 ||p|| w e9 +^(l|p|||,llv / p7M^i|li 2 



"* 1 37^£ 1 _ 2(39-67! 

-Hp « 

+ ||p|iin V7(4^ • v^iiH, + hvTWp" • v(a)- 1 v(p(p))||| p + iiV/Mp^iiL 



+ [lv7Wpft*), 

where we have 9 which depends of t for the bootstrap. 
Estimate on the term J * f TN (pg — ^pd s v) f ' (s)d s v\dxds 

For the term J * fj N (pg — ^pd s v)f(s)d s v\dxds, we proceed exactly as in the previous 
section. 

Control on the norm p G L°° 

Next we come back to equation (|3.36j) to get a control in norm L°° on the density, more 
precisely we have: 

log(p(i,x)) <log(||po||L-) + C||(A)- 1 divm ||L- +C||(A)- 1 div(pu)|| L ~ 

ft 1 (7.88) 
+ C / —=\\[uj, RiRj](pui)(s,-)\\ L °ods. 



We recall then from the previous section that f(s) Av± G L 2 (L P ) with | = | + by 
Gagliardo-Nirenberg inequality we have: 

||V^||1 91 < \\^f{s)Av 1 f^ ) \\Vv 1 \\% 

with ^" = ^ + ^ — — ^ ~ Tv an d 9 > Y' ^ e nave then by using the results of R. 
Coifman et al in |19j : 

(m)^[(v 1 ) j ,R i R j ](pu i )(8,')\\ w x, a < \\V Vl \\l4yW)^n\\P\\p^n\\ L:p ^_ 
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q 

1 ., 1 2 . 



x IIpII v , 



(7.89) 

with i = 22^2 + 37^2+1+^ + 1 < 1 let 7g 2 _ 3g(5 + + lg7 > Q After a 
small calculus we obtain <? > jf + yj;7 + ^y / 97 2 — 267 + 23. From (|7.89p . we have 
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(/( S ))t-?||[(t;i) i ,/J i i2 7 -](pu i )( S ,.)[| W i,a G L 2 . 
We obtain then by Sobolev embedding: 

, 1 - \ J (s) 1 *-H[(v 1 ) 3 ,R l R j ]{pu i ){s,-)\\ L ~ds\ <C\\p^u\\ Toc .« , 2 • 

We proceed similarly for the term \\[(A~ 1 (P(p))j, RiRj](pui)(s, ■) ||z,i(z,°°) - This term is 
crucial because it gives the value of q that we must choose. Indeed we have the results 
of R. Coifman et al in |19| if q < 3j: 

\\[(A-\P(p)) j ,R i R j ](pu i )(8,')\\ W i,i> < \\P( P )Li\\\\p^u\\ «\\p\\ q 2 , 

with ^ = + 3 7~Q _|_ 2 < I which implies 3q 2 — q(9*y + 5) + 67 > 0. However we see 

by the previous inequality that we must have q > 37. 

1 

As uq G L°°, we can show that ppu G L°°(L P ) with p arbitrarly large because A = 0. We 

1 , 1 -fx 

have then ppu G L°°(L P ), p p G L°°(L p ). We conclude that by Holder's inequalities 

that: 

g 1 1 

- + - < 

7^3 

let 5 > 3(7 + 1). 

By this way we control log p G L°° by the fact that q > 3(7 + 1) and that (|7.85p and 
(|7.87p are verified. It means that p G L°° and = G L°°, from theorem 11.21 we have seen 
that the inequality (|1.10p and (jl.lip are preserved during the time if p G L°° . We can 
now assume that (p, u) verify (jl.lOp and (jl.lip . We want now prove the uniqueness of 
this solution. 

7.2 Uniqueness 

You want prove now the result of uniqueness. To do it we want use the result of P. 
Germain in [3T] . In the sequel we will note (p\ , u\ ) the solution of the theorem 11.21 which 
exits on the time interval [0,To]. We have shown that our solution check p G L 00 (L 00 ). 
By theorem 11.21 we obtain that our solution verify the following inequalities: 

sup / ip(t,x)\u(t,x)\ 2 + \P(p(t,x))\ + a(t)\Vu(t,x)\ 2 dx 

0<t<+oo JjN 2 

+ sup / [\p{t,x)f{t) N {p\u{t,x)\ 2 + \Vu{t,x)\ 2 )dx 

0<t<+ooJTN I (7.90) 

r+00 r- 

+ / / [\Vu\ 2 + f{s)p\u\ 2 + \uj\ 2 ) + a N \Vii\ 2 ]dxdt 
Jo Jt n 

< C(C + C f )\ 

and we obtain moreover: 

{ ^d t ueL 2 (L 2 (T N )), 
VtVueL 2 T (H 2 (T N )), 

VtG = Vt[{\ + 2/i)divu — P(p)] G L|.(i/ 1 (T Ar )), l ' ' 

k VtVu G (L 2 (T JV )), 
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Now by the result of P. Germain in [31], we are able to prove that (p,u) = (p\,u\) 
on [0,To]. To see this we have just to verify that (p\,u\ verify the conditions of the 
theorem 2.2 of [31]. For simplicity we prove only the result for N = 3. We know that 
Vpi G L°°(B e Nl ) ^ L^(L N ) and Vu t G L l To {L°°). The main thing is to see that 
\fiu\ G Lj, Q (L N ). 

We recall by Gagliardo-Nirenberg inequalities that: 

V*||wi||i3 < (t' _ ^||«i||i,2)3(t5~f||Vtti||£a)5tf . 
From the inequalities (|1.14|) . we deduce that \ftii\ G L 2 (L 3 ). 

We have proved that our solution is unique on [0, To], we have to see now what happen 
for t > Tq. When t > Tq, we have in this case by using inequalities (|7.90[> . 

sup / [~p(t,x)\u(t,x)\ 2 + \P(p(t,x))\ + a(t)\Vu(t,x)\ 2 dx 
T <t<+ooJr N * 

+ sup / [\p{t,x)f{T Q ) 2 {p\ii(t,x)\ 2 + \Vu{t,x)\ 2 )dx 

T <t<+ooJjN 2 (7.92) 

r+oc r 

+ / / [|Vn| 2 + /(T )) /0 |^| 2 + k| 2 ) + /(To) 2 |Vn| 2 ](ixdt 

<C(Co + C f f, 

We can check easily that v\ G T^ oc ([To, +oo), B^^) so by proceeding similarly as sub- 
subsection ESM we obtain that p G L°° [R, B 1 ^) . We get then that Vu G L 1 (E, L°°) 
and we conclude by using again the result of P. Germain in [31J. 

References 

[1] H. ABIDI and M. PAICU. Fquation de Navier-Stokes avec densite et viscosite vari- 
ables dans l'espace critique. Annates de Vinstitut Fourier, 57 no. 3 (2007), p. 883-917. 

[2] H. BAHOURI and J.-Y. CHEMIN, Equations d'ondes quasilineaires et estimation de 
Strichartz, Amer. J. Mathematics, 121, (1999), 1337-1377. 

[3] H. BAHOURI, J.-Y. CHEMIN and R. DANCHIN. Fourier analysis and nonlinear 
partial differential equations, to appear in Springer. 

[4] J.-M. BONY, Calcul symbolique et propagation des singularites pour les equations 
aux derivees partielles non lineaires, Annales Scientifiques de Vecole Normale 
Superieure, 14, (1981), 209-246. 

[5] D. BRESCH and B. DESJARDINS, On the construction of approximate solutions 
for the 2d viscous shallow water model and for compressible Navier-Stokes models. J. 
Maths Pures et AppL, (86), 4, (2006), 262-268. 

[6] D. BRESCH and B. DESJARDINS, Existence of global weak solutions for a 2D Vis- 
cous shallow water equations and convergence to the quasi-geostrophic model. Comm. 
Math. Phys., 238 (1-2), 211-223, 2003. 



34 



[7] D. BRESCH and B. DESJARDINS, Existence of global weak solutions to the Navier- 
Stokes equations for viscous compressible and heat conducting fluids, Journal de 
Mathematiques Pures et Appliques, Volume 87, Issue 1, January 2007, Pages 57-90. 

[8] D. BRESCH and B. DESJARDINS, Some diffusive capillary models of Kortweg type. 
C. R. Math. Acad. Sci. Paris, Section Mecanique, 332(11), 881-886, 2004. 

[9] D. BRESCH, B. DESJARDINS and C.-K. LIN, On some compressible fluid mod- 
els: Korteweg, lubrification and shallow water systems. Comm. Partial Differential 
Equations, 28(3-4), 843-868, 2003. 

[10] L. CAFFARELLI, R. KOHN and L. NIRENBERG, Partial regularity of suitable 
weak solutions of the Navier-Stokes equations., Comm. Pure Appl. Math., 35(6): 771- 
831, 1982. 

[11] F. CHARVE and R. DANCHIN, Global existence in critical spaces for compressible 
Navier-Stokes equations, preprint. 

[12] J.-Y. CHEMIN, Theoremes d'unicite pour le systeme de Navier-Stokes tridimension- 
nel, J. d Analyse Math. 77, (1999), 27-50. 

[13] J.-Y. CHEMIN and N. LERNER, Flot de champs de vecteurs non lipschitziens et 
equations de Navier-Stokes, J. Differential Equations, 121, (1992), 314-328. 

[14] Q. CHEN, C. MIAO and Z. ZHANG, Global well-posedness for the compressible 
Navier-Stokes equations with the highly oscillating initial velocity, arXiv:0907.4540v2. 

[15] Q. CHEN, C. MIAO and Z. ZHANG, Well-posedness in critical spaces for 
the compressible Navier-Stokes equations with density dependent viscosities, arXiv 
:0811.4215vl November 2008. 

[16] Y. CHO and H. KIM, Unique solvability for the density-dependent Navier-Stokes 
equations. 

[17] H. J. CHOE and H. KIM, Strong solution of the Navier-Stokes equations for isen- 
tropic compressible fluids, J. Differential Equations, 190, (2003), 504-523. 

[18] H. J. CHOE and H. KIM, Strong solution of the Navier-Stokes equations for non- 
homogeneous incompressible fluids, Math. Meth. Appl. Sci. 28, (2005), 1-28. 

[19] R. COIFMAN, P.-L. LIONS, Y. MEYER and S. SEMMES, Compensated- 
compactness and Hardy spaces, J. Math. Pures. Appl, 72 (1993), p247-286. 

[20] F. COQUEL, D. DIEHL, C. MERKLE and C. ROHDE, Sharp and diffuse interface 
methods for phase transition problems in liquid-vapour flows. Numerical Methods for 
Hyperbolic and Kinetic Problems, 239-270, IRMA Lect. Math. Theor. Phys, 7, Eur. 
Math. Soc, Zurich, 2005. 

[21] R. DANCHIN, Fourier analysis method for PDE's, Preprint, Novembre 2005. 



35 



[22] R. DANCHIN, Local Theory in critical Spaces for Compressible Viscous and Heat- 
Conductive Gases, Communication in Partial Differential Equations, 26 (78), 1183- 
1233, (2001). 

[23] R. DANCHIN, Global Existence in Critical Spaces for Flows of Compressible Viscous 
and Heat-Conductive Gases, Arch. Rational Mech.Anal.160, (2001), 1-39. 

[24] R. DANCHIN, On the uniqueness in critical spaces for compressible Navier-Stokes 
equations. NoDEA Nonlinear Differentiel Equations Appl, 12(1):111-128, 2005. 

[25] R. DANCHIN, Well-Posedness in Critical Spaces for Barotropic Viscous Flu- 
ids with Truly Not Constant Density, Communications in Partial Differential 
Equations,32:9, 1373-1397. 

[26] B. DESJARDINS, Regularity of weak solutions of the compressible isentropic 
Navier-Stokes equations, Comm Partial Differential Equations 22 (1997) 977-1008. 

[27] E. FEIREISL, Dynmamics of Viscous Compressible Fluids, Oxford Lecture Series in 
Mathematics and its Applications, 26. 

[28] E. FEIREISL. Compressible Navier-Stokes equations with a non-monotone pressure 
law. J. Differential Equations, 184(1), 97-108, 2002. 

[29] E. FEIREISL. On the motion of a viscous, compressible, and heat conducting equa- 
tion. Indiana Univ. Math. J., 53(6), 1705-1738, 2004. 

[30] E. FEIREISL, A. NOVOTNY and H. PETZELTOVA. On the existence of globally 
defined weak solutions to the Navier-Stokes equations. J. Math. Fluid Mech., 3(4), 
358-392, 2001. 

[31] P. GERMAIN, Weak stong uniqueness for the isentropic compressible Navier-Stokes 
system, preprint. 

[32] Z. GUO, Q. JIU and Z. XIN. Spherically symmetric isentropic compressible flows 
with density-dependent viscosity coefficients. SI AM J. Math. Anal. (2008), 1402-1427. 

[33] B. HASPOT, Cauchy problem for viscous shallow-water equations with a term of 
capillarity, accepted in MS AS. 

[34] B. HASPOT, Cauchy problem for viscous shallow water equations with a term of cap- 
illarity , accepted in Hyperbolic Problems: Theory, Numerics and Application PSAPM 
volume, edited by Drs. Tadmor, Liu, Tzavaras. 

[35] B. HASPOT, Local well-posedness results for density-dependent incompressible flu- 
ids, i4rw«,:0902.1982 (February 2009). 

[36] B. HASPOT, Well-posedness in critical spaces for barotropic viscous fluids , 
Areiv,:0902.1982 (February 2009). 

[37] D. HOFF. Uniqueness of weak solutions of the Navier-Stokes equations of multidi- 
mensional compressible flows. SI AM. J. Anal.,vol 37, No6, 1742-1760 (2006). 



36 



[38] D. HOFF. Compressible flow in a Half-space with Navier boundary conditions, Jour- 
nal of Mathematical Fluid Mechanics, Volume 7, Number 3 / August, (2005) 315-338. 

[39] D. HOFF. Dynamics of singularity surfaces for compressible, viscous flows in two 
space dimensions, Communications on Pure and Applied Mathematics, Volume 55 Issue 
11, 1365 - 1407 (2002). 

[40] D. HOFF. Global existence for ID, compressible, isentropic Navier-Stokes equations 
with large initial data. Trans. Amer. Math. Soc, 303(1), 169-181, 1987. 

[41] D. HOFF. Global well-posedness of the Cauchy problem for nonisentropic gas dy- 
namics with discontinuous data, J. Diff. Eq., 95 (1992), p33-73 

[42] D. HOFF. Discontinuous solutions of the Navier-Stokes equations for multidimen- 
sional flows of the heat conducting fluids. Arch. Rational Mech. Anal., 139, (1997), p. 
303-354. 

[43] D. HOFF. Global solutions of the Navier-Stokes equations for multidimensional 
compressible flow with discontinuous initial data. J. Differential Equations, 120(1), 
215-254, 1995. 

[44] D. HOFF. Strong convergence to global solutions for multidimensional flows of com- 
pressible, viscous fluids with polytropic equations of state and discontinuous initial 
data. Arch. Rational Mech. Anal, 132(1), 1-14, 1995. 

[45] D. HOFF and K. ZUMBRUM. Multi-dimensional diffusion waves for the Navier- 
Stokes equations of compressible flow, Indiana University Mathematics Journal, 1995, 
44, 603-676. 

[46] S. JIANG and P. ZHANG. Axisymetrics solutions of the 3D Navier-Stokes equations 
for compressible isentropic fluids. J. Math. Pures Appl. (9), 82(8): 949-973, 2003. 

[47] A. V. KAZHIKOV. The equation of potential flows of a compressible viscous fluid for 
small Reynolds numbers: existence, uniqueness and stabilization of solutions. Sibirsk. 
Mat. Zh., 34 (1993), no. 3, p. 70-80. 

[48] A. V. KAZHIKOV and V. V. SHELUKHIN. Unique global solution with respect to 
time of initial-boundary value problems for one- dimensional equations of a viscous 
gas. Prikl. Mat. Meh., 41(2): 282-291, 1977. 

[49] H.-L. LI, J. LI and Z. XIN. Vanishing of vacuum states and blowup phenomena of 
the compressible Navier-Stokes equation. Accepted in Comm. Math. Phys, (2007). 

[50] P.-L. LIONS. Mathematical Topics in Fluid Mechanics, Vol 2, Compressible models, 
Oxford University Press, (1998). 

[51] A. MATSUMURA and T. NISHIDA. The initial value problem for the equations of 
motion of compressible viscous and heat-conductive fluids. Proc. Japan Acad. Ser. A 
Math. Sci, 55(9), 337-342, 1979. 



37 



[52] A. MELLET and A. VASSEUR, On the isentropic compressible Navier-Stokes equa- 
tion, Arxiv preprint math.AP/0511210, 2005 - arxiv.org 

[53] A. MELLET and A. VASSEUR, L p estimates for quantities advected by a com- 
pressible flow, Journal of Mathematical Analysis and Applications, Volume 355, Issue 
2, 548-563 (2009). 

[54] A. MELLET and A. VASSEUR, Existence and uniqueness of global strong solutions 
for one-dimensional compressible Navier-Stokes equations. SI AM J. Math. Anal. 39 
(2007/2008), no. 4, 1344-1365. 

[55] J. NASH, Le probleme de Cauchy pour les equations differentielles d'un fluide 
general, Bulletin de la Societe Mathematique de France, 1962, 90, 487-497. 

[56] A. NOVOTNYY and I. STRASKRABA, Introduction to the mathematical theory 
of compressible flow, Oxford lecture series in mathematics and its applications. 27. 

[57] C. ROHDE, On local and non-local Navier-Stokes-Korteweg systems for liquid- 
vapour phase transitions. ZAMMZ. Angew. Math. Mech., 85(2005), no. 12, 839-857. 

[58] T. RUNST and W. SICKEL. Sobolev spaces of fractional order, Nemytskij operators, 
and nonlinear partial differential equations, de Gruyter Series in Nonlinear Analysis 
and Applications, 3. Walter de Gruyter and Co., Berlin (1996). 

[59] R. SALVI and I. STRASKRABA, Global existence for viscous compressible fluids 
and their behavior at t — > +oo, J. Fac. Sci. Univ. Tokyo Sect.lA, Math, 40 (1993) 
17-51. 

[60] D. SERRE. Solutions faibles globales des equations de Navier-Stokes pour un fluide 
compressible., 303(13): 639-642, 1986 

[61] V. A. SOLONNIKOV. Estimates for solutions of nonstationary Navier-Stokes sys- 
tems. Zap. Nauchn. Sem. LOMI, 38, (1973), p.153-231; J. Soviet Math. 8, (1977), p. 
467-529. 

[62] V. A. SOLONNIKOV. Solvability of the initial boundary value problem for the 
equation of a viscous compressible fluid, J. Sov. Math., 14 (1980), pll20-1133. 

[63] V. VALLI, W. ZAJACZKOWSKI. Navier-Stokes equations for compressible fluids: 
global existence and qualitative properties of the solutions in the general case. Com- 
mun. Math. Phys., 103, (1986) no 2, p. 259-296. 

[64] J. D. VAN DER WAALS. Thermodynamische Theorie der Kapillaritat unter Vo- 
raussetzung stetiger Dichteanderung, Z. Phys. Chem. 13, (1894), 657-725. 

[65] A. VASSEUR. A new proof of partial regularity of solutions to Navier-Stokes equa- 
tions. NoDEA Nonlinear Differential Equations Appl, 14(5-6):753-785, 2007. 

[66] V. A. WEIGANT, Example of non-existence in the large for the problem of the 
existence of solutions of Navier-Stokes equations for compressible viscous barotropic 
fluids, Dokl. Akad. Nauk 339 (1994) 155-156 (in Russian). 



38 



[67] W. WANG and C-J XU. The Cauchy problem for viscous shallow water equations. 
Rev. Mat. Iberoamericana 21., no. 1 (2005), 1-24. 



39 



